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ABSTRACT 

Standard shot-noise models, which seek to explain the broadband noise variability 
that characterises the X-ray lightcurves of X-ray binaries and active galaxies, predict 
that the power spectrum of the X-ray lightcurve is stationary (i.e. constant amplitude 
and shape) on short time-scales. We show that the broadband noise power spectra of 
the black hole candidate Cyg X-l and the accreting millisecond pulsar SAX J1808.4- 
3658 arc intrinsically non-stationary, in that RMS variability scales linearly with flux. 
Flux-selected power spectra confirm that this effect is due to changes in power-spectral 
amplitude and not shape. The lightcurves of three Seyfert galaxies are also consistent 
with a linear relationship between RMS variability and flux, suggesting that it is an 
intrinsic feature of the broadband noise variability in compact accreting systems over 
more than 6 decades of central object mass. The RMS variability responds to flux 
variations on all measured time-scales, raising fundamental difficulties for shot-noise 
models which seek to explain this result by invoking variations in the shot parameters. 
We suggest that models should be explored where the longest time-scale variations 
are fundamental and precede the variations on shorter time-scales. Possible models 
which can explain the linear RMS-flux relation include the fractal break-up of large 
coronal flares, or the propagation of fluctuations in mass accretion rate through the 
accretion disk. The linear relationship between RMS variability and flux in Cyg X-f 
and SAX Jf808. 4-3658 is offset on the flux axis, suggesting the presence of a second, 
constant-flux component to the lightcurve which contributes ~ 25% of the total flux. 
The spectrum of this constant component is similar to the total spectrum, suggesting 
that it may correspond to quiet, non- varying regions in the X-ray emitting corona. 

Key words: X-rays: stars - stars: individual: Cyg X-f - stars: individual: 
SAX Jf 808. 4-3658 - galaxies: active - galaxies: Seyfert 



1 INTRODUCTION 

The broadband noise variability which characterises the X- 
ray lightcurves of both X-ray binary systems (XRBs) and 
active galactic nuclei (AGN) is commonly modelled as a 
shot-noise process. The simplest shot-noise models, where 
the lightcurve is constructed from a stochastic series of inde- 
pendent overlapping shots with a single characteristic decay 
time-scale (e.g. Terrell f 972) , produce a red-noise type power 
spectrum at high frequencies (P(u) oc v~ a , where a = 2), 
which flattens to a — at low frequencies (i.e. the variabil- 
ity becomes white-noise on long time-scales). This simple 
form is qualitatively similar to the power-spectral shapes 
of XRBs, which flatten below ~ O.f Hz (e.g. van der Klis 
1995) and AGN, which flatten on much longer time-scales 
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(< 10" 6 Hz, M c Hardy 1988, Edelson & Nandra 1999), but 
fails to reproduce the typically observed red-noise slopes, 
a ~ 1.0. More complex shot-noise models, which invoke 
a broad distribution of decay time-scales, can reproduce 
the correct high-frequency shapes (Lehto 1989; Lochner, 
Swank & Szymkowiak 1991), including the additional high- 
frequency break which is seen above ~ 1 Hz in some XRB 
power spectra (e.g. Cyg X-l, Belloni & Hasinger 1990). 
Physical interpretations of shots typically invoke magnetic 
flares in the X-ray emitting corona which is likely to be 
common to both XRBs and AGN (e.g. Poutanen & Fabian 
1999). 

It is well-known that the power spectra of XRBs are non- 
stationary (i.e. have a variable shape and amplitude) on long 
time-scales (>days, e.g. Belloni & Hasinger 1990). There is 
some evidence for non-stationarity of AGN power spectra 
on much longer time-scales (Uttley et al. 1999). The shot- 
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Figure 1. Local flux dependence of mean a for Cyg X-l and SAX J1808. 4-3658. The dotted lines mark the best-fitting linear models 
described in the text. 



noise models described above can account for non-stationary 
power-spectra by invoking changes in the parameters of the 
shot-noise process on long time-scales, but predict that on 
short time-scales the power spectrum should be stationary. 
Because individual power-spectra measured on very short 
time-scales (seconds) are noisy, this prediction is difficult 
to test. However, we can test whether power-spectral shape 
and/or amplitude varies systematically with short time-scale 
variations in X-ray flux by binning the power-spectra of 
small segments of observed XRB and AGN lightcurves ac- 
cording to their flux. 

The relationship between long-time-scale variations in flux 
and power-spectral amplitude has been studied in a few 
XRBs, e.g. a somewhat complex correlation between RMS 
variability and flux has been noted during the decay of the 
hard X-ray transient GRO J0422+32 (Denis et al. 1994). 
Surprisingly, the possibility of systematic short-time-scale 
variations of power-spectral shape and amplitude with flux 
has not previously been investigated in detail. The stan- 
dard practice of normalising power spectra by squared mean 
flux (and RMS variability by mean flux), for comparison 
between different objects and instruments, causes informa- 
tion regarding the flux to be lost. In this letter, we show 
that this practice has helped to keep hidden a startling fact 
about the X-ray variability of XRBs and AGN: that the 
RMS variability of their broadband noise component scales 
with flux variations on time-scales as short as seconds, while 
the power-spectral shape remains constant. This discovery 
poses particular problems for standard shot-noise models, as 
we shall discuss. 



2 THE FLUX-DEPENDENT VARIABILITY OF 
X-RAY BINARIES 

We obtained public archival data from Rossi X-ray Timing 
Explorer (RXTE) observations of Cyg X-l (in the low state), 
observed 23 October 1996 (Nowak et al. 1999) and the re- 
cently discovered accreting millisecond pulsar, SAX J1808.4- 



3658 (Wijnands & van der Klis 1998a), observed 18 April 
1998, of duration 18 ksec and 23 ksec respectively (using 
the PCA instrument and including only periods where all 
5 Proportional Counter Units (PCUs) were switched on). 
Cyg X-l is an obvious choice for study, as it is the most 
well-known black hole X-ray binary. The accreting millisec- 
ond pulsar SAX J1808. 4-3658 is a good example of a neutron 
star system which displays strong broadband noise variabil- 
ity, so any similarity in the variations of the power spectra 
of these objects will provide strong evidence that the power- 
spectral variability is intrinsic to the process that produces 
the broadband noise, independent of whether the central ob- 
ject has a hard surface or an event horizon. Furthermore, the 
observations we select show no evidence for discernable low- 
frequency quasi-periodic oscillations (QPOs) or components 
other than the broadband noise which we wish to investigate 
here. Using PCA binned and event mode data from the ob- 
servations of Cyg X-l and SAX J1808.4-3658 respectively, 
we made lightcurves with a resolution of 16 ms for each 
source in the 2-13.1 keV energy range. The mean fluxes in 
this energy range are 3377 count s _1 and 409 count s _1 for 
Cyg X-l and SAX J1808.4-3658 respectively. 
We first investigated the dependence of RMS variability on 
local X-ray flux, in other words, does the RMS variability of 
a small section of lightcurve depend on the mean flux of that 
section of lightcurve? In order to answer this question, wc 
split each lightcurve into 10 s segments and measured the 
power spectrum for each segment (subtracting the contri- 
bution due to Poisson noise estimated from the mean total 
count rate of the lightcurve segment). We applied the stan- 
dard RMS-squared normalisation to the power spectrum, 
so that integrating the power over a given frequency range 
yields the contribution to lightcurve variance due to varia- 
tions in that frequency range. Taking the square root of the 
integrated power yields the RMS variability a. By using this 
method, rather than simply measuring a directly from the 
lightcurve, we can study how the amplitude of the power 
spectrum in a defined frequency range responds to changes 
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Figure 2. Flux-dependent power spectra of Cyg X-l, correspond- 
ing to the first and fourth flux quartiles (see text for details). 



in flux. We first chose to measure a over the frequency range 
0.1-10 Hz, which incorporates the full spectrum of broad- 
band noise from white-noise to red-noise, as can be seen 
from the power spectra of these observations (Nowak et al. 
1999, Wijnands & van der Klis 1998b). Due to the stochas- 
tic nature of red-noise lightcurves, there is significant scatter 
in individual measures of a, so we binned a as a function 
of segment flux using relatively narrow flux bins to exam- 
ine the form of the <r-flux relationship over a broad range 
of flux. Using a minimum number of 30 measurements per 
bin to ensure an accurate estimate of the standard error in 
the mean a, the resulting er-flux relation for both sources 
is plotted in Fig. ^. The lightcurves of both Cyg X-l and 
SAX J1808. 4-3658 show a remarkably linear dependence of 
RMS variability on flux. 

Visual inspection of Fig. [j] shows that the er-flux trend in 
both sources does not pass through the origin, i.e. there 
is a constant offset in flux and possibly also a. To test 
the goodness of fit of a linear model, including a constant 
offset from the origin, we fitted a function of the form 
a — k(F—C), where F is the flux and k and C are constants. 
This linear model provides a good fit to the data, yielding 
X 2 values of 28.7 (for 22 degrees of freedom) for Cyg X- 
1 and 17.2 (20 degrees of freedom) for SAX J1808.4-3658. 
The best fitting model parameters were k = 0.326 ± 0.017, 
C = 850 ±130 count s" 1 for Cyg X-l and k = 0.305 ±0.026, 
C = 81 ± 30 count s _1 for SAX J1808.4-3658 (uncertainties 
are at the 90% confidence level for 2 interesting parameters) . 
Note that the gradient of the a-F trend k, is equivalent to 
the fractional RMS variability of the variable-cr component 
of the lightcurve. The constant C represents a second com- 
ponent to the lightcurve, which does not follow the linear 
a-F trend but may still contribute to the total value of a. It 
is not possible to disentangle the mean flux level of this sec- 
ond, constant-cr component and the contribution it makes 
to the total a, however the value of C represents the flux of 
this component in the limit where it does not vary. 

Changes in RMS variability measured over a given fre- 
quency range may be associated with changes in power 
spectral amplitude and/or shape. For example, Belloni & 
Hasinger (1990) show that the fractional RMS variability 
of Cygnus X-l is well correlated with changes in the po- 



sition of the low-frequency break in the power spectrum, 
but note that these variations occur on long time-scales and 
show no relationship to X-ray flux. It seems unlikely that 
the strict linear relationship we have found between RMS 
variability and flux could be caused by changes in power- 
spectral shape over the entire 0.1-10 Hz frequency range, 
because the power spectrum over this range is not a simple 
power law. Instead, the simplest explanation is that only 
the power-spectral amplitude changes. We can confirm this 
interpretation of the cr-flux relation by plotting power spec- 
tra which are binned according to lightcurve flux. Fig. [| 
shows flux-binned power spectra of Cygnus X-l, obtained 
by averaging the power spectra of lightcurve segments of 
32 s duration, according to whether the segment flux occu- 
pies the first or fourth quartile of the overall distribution of 
fluxes. The resulting power spectra correspond to fluxes of 
3092 count s _1 and 3679 count s _1 for the first and fourth 
flux quartiles respectively, and show that the amplitude of 
the broadband noise variability responds to flux in the same 
way across the entire power spectrum. There does not ap- 
pear to be any systematic change in power spectral shape. 
We obtain a similar result for SAX J1808.4-3658. 



3 THE FLUX-DEPENDENT VARIABILITY OF 

AGN 

Given the strong similarities between the power-spectral 
shapes of AGN and XRBs, which suggest that the same 
variability processes are at work, independent of the mass of 
the central object, we might expect AGN to show the same 
flux-dependent X-ray variability as XRBs. To test this pos- 
sibility, we obtained quasi-continuous lightcurves of ~days 
duration for three Seyfert galaxies, NGC 4051, NGC 5506 
and MCG-6-30-15, observed in Dec 1996, Jun 1997 and Aug 
1997 respectively and available in the RXTE public archive. 
Because PCUs 3 and 4 were switched off for much of the 
duration of each observation, we only use data from PCUs 
0, 1 and 2. We make the lightcurves from Standard 2 data, 
in the 2-10 keV energy band, using data from only the top 
layer of each PCU in order to minimise the contribution of 
instrumental background which adds to the Poisson noise 
level, applying standard selection criteria (e.g. as described 
in Uttley 1999) to extract good time intervals, and estimat- 
ing background lightcurves using the L7 model. The obser- 
vations are of total useful exposures ~ 75 ks, 90 ks and 
340 ks for NGC 4051, NGC 5506 and MCG-6-30-15 respec- 
tively. We made power spectra from continuous lightcurve 
segments of length ~ 2500 s (corresponding to periods be- 
tween Earth occultations of the source), and measured the 
Poisson-noise-subtracted integrated power in the 5 x 10 -3 - 
5 x 10 -4 Hz frequency range. 

Due to the faint nature and lower variability (relative to 
XRBs) of AGN on the time-scales sampled, we average the 
integrated power into two flux bins, corresponding to mean 
segment fluxes lying below or above the mean flux of the en- 
tire observation. The resulting integrated powers (including 
standard errors) are shown in Table |l|. Note that we do not 
directly determine the RMS variability, a for each segment, 
because statistical fluctuations in the true Poisson noise 
level, combined with the relatively small amount of high- 
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Table 1. Flux-related changes in the variance of NGC 4051, NGC 5506 and MCG-6-30-15. 

Low flux High flux 
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NGC 4051 


3.4 


1.9 


17 


0.047 ±0.024 


11.4 % 


5.3 


16 


0.32 ± 0.07 


10.7 % 


NGC 5506 


28.1 


25.5 


20 


0.21 ±0.05 


1.8 % 


31.2 


19 


0.41 ± 0.11 


2.1 % 


MCG-6-30-15 


12.2 


10.3 


73 


0.28 ± 0.05 


5.1 % 


14.1 


70 


0.56 ±0.07 


5.3 % 



Mean data are shown for low and high flux segments (corresponding to segment fluxes below and above the mean for the entire 
observation, given by fitot)- fJ- is the mean flux (2-10 keV, count s — 1 ) for each flux bin, n is the number of segments in the flux bin, a 1 
is the mean variance (count 2 s -2 ) and o-f rac is the fractional RMS variability (o"f rac = ( a ' 2 /A t2 ) 5 ) 



frequency power seen in AGN lightcurves, sometimes lead 
to negative integrated powers after subtraction of the noise 
estimate. Segments with negative integrated power cannot 
be used to determine a, but can be used to determine the 
mean a 2 and its error (note that the fractional RMS can be 
estimated from the mean a 2 and is also shown in Table |l|) . 
These data clearly show that the variance of the X-ray 
lightcurves of all three AGN is dependent on X-ray flux, 
and furthermore the relationship between the RMS vari- 
ability and the mean segment flux seems to be linear, such 
that the fractional RMS variability remains approximately 
constant, despite significant flux changes (e.g. ~ factor 3 for 
NGC 4051). There is not sufficient data to confirm the exis- 
tence of a constant flux component to the AGN lightcurves, 
of the same relative strength as that seen in the XRBs, al- 
though we note that the existence of such a component is 
not ruled out. 
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Figure 3. Energy-dependent flux of possible constant compo- 
nents to the lightcurves of Cyg X-l and SAX J1808.4-3658, rela- 
tive to the mean flux in each band. 



4 DISCUSSION 

The relation between flux and the amplitude of broadband 
noise variability which we have presented here, suggests that 
the lightcurves of XRBs (and possibly AGN) are made from 
at least two components. One component shows a strik- 
ing linear dependence of RMS variability on flux, while the 
other component may contribute a constant RMS to the 
lightcurve or, more simply, may not vary at all. We can 
crudely estimate the spectral shape of this second, possi- 
bly constant component by determining the RMS-flux rela- 
tion for lightcurves in four energy bands and measuring the 
value of the intercept on the flux axis (C) for each band. In 
Fig. ^ we plot the value of C in each energy band expressed 
as a fraction of the mean flux in that band. The errors on 
the relative flux are determined from errors in C estimated 
from fitting a linear model to the RMS vs. flux data for 
that energy band. Note however that because the fluxes in 
each energy band are strongly correlated with one another 
(e.g. Maccarone, Coppi & Poutanen 2000), the relative un- 
certainty between energy bands is likely to be significantly 
smaller than indicated by the error bars we determine here. 
With this caveat in mind, Fig. ^ suggests that, if C repre- 
sents the flux of a constant component to the lightcurve, its 
spectral shape is very similar to that of the total spectrum 
and hence the variable component. 

This result seems strange if we naively expect that any con- 
stant components to the lightcurve should have a different 
origin and therefore a different spectrum to the variable com- 
ponent. We note however, that provided that the tempera- 
ture of the corona is maintained, it will continue to Comp- 
tonise seed photons regardless of whether it is dynamically 



variable (e.g. flaring) or 'quiet' and non-variable. Therefore 
the constant component to the X-ray lightcurve may be 
associated with quiet regions of the corona or patches of 
the corona above inactive regions of the accretion disk (if 
variability is driven by variations in seed photon number). 
The ~ 25% of total flux which the constant component con- 
tributes could indicate a ~ 25% covering fraction of quiet 
coronal regions (this may represent a time-averaged cover- 
ing fraction, since coronal regions might switch from being 
quiet to being variable) . An alternative possibility is that the 
constant component is due to Thomson scattering of the pri- 
mary continuum photons in extended, hot gas with radius 
R > 10 light-seconds to smear out any variability, tempera- 
ture T ~ 10 keV so that no absorption edges are observed 
and optical depth r ~ 0.25. However, a plasma of this extent 
should be a strong source of thermal bremsstrahlung emis- 
sion. Assuming a spherical geometry and uniform density, 
the bremsstrahlung luminosity from the extended plasma 
region would scale linearly with its radius so that such a 
region would contribute > 10% of the X-ray luminosity of 
the system in bremsstrahlung emission and is ruled out by 
observations. 

The variable component of the lightcurve shows a linear 
dependence of RMS variability on flux, which is not con- 
sistent with the expectations of simple shot-noise models, 
which predict that the power spectrum is stationary. At first 
glance, the obvious explanation is that the parameters of the 
shot-noise process are changing on short time-scales. For ex- 
ample, if the average amplitude of shots varies on short time- 
scales, the flux and RMS variability will vary in proportion 
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duced on time-scales greater than the viscous time-scale, so 
that longer time-scale variations are produced in the outer 
disk and shorter time-scale variations in accretion rate are 
superimposed on them as they propagate inwards. Thus, a 
linear RMS-flux relation will be produced if the fractional 
amplitude of accretion rate variations is independent of the 
actual accretion rate (as suggested by the model, which as- 
sumes that changes in accretion rate are caused by fractional 
fluctuations in the disk viscosity parameter) . Either of these 
suggested models might, in principle, explain the RMS-flux 
relation we see, but we leave a detailed comparison of their 
predictions with observations for a future work. 
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Figure 4. Flux-dcpcndcncc of mean a (2-20 Hz band) for Cyg X- 
1, determined from 1 s lightcurve segments. 



to the shot amplitude. Hence the linear relation between 
RMS variability and flux reflects the linear relationship of 
both variables to a single underlying varying parameter, the 
shot amplitude. However, this model cannot simply explain 
why the RMS variability is dependent on flux on all mea- 
sured time-scales, a fact which becomes apparent if we plot 
the RMS-flux relation for Cygnus X-l using 1 s time bins, de- 
termining o in the 2-20 Hz frequency range (see Fig. ^J). The 
fact that the linear relationship between a and flux extends 
to a greater range of fluxes than are obtained by determin- 
ing the relation using 10 s segments, implies that the linear 
RMS-flux relation applies to flux variations on time-scales 
as low as 1 s (which are averaged out by the 10 s segments 
used earlier). Therefore, if shot parameter variations are the 
cause of the RMS-flux relation, they must occur on as broad 
a range of time-scales as the variations in the lightcurve it- 
self (at least down to 1 s), and we are left with the circular 
problem of trying to explain the broadband noise in X-ray 
lightcurves with a model which itself requires a broadband 
noise component (to describe the variation in shot parame- 
ters) in order to be consistent with the data. 
The simplest solution may be to discard the standard shot- 
noise models altogether. Instead, the linear dependence of 
RMS variability on flux may be more simply explained if the 
lightcurves are instead made 'from the top down', with the 
primary cause of variability being large, long-time-scale vari- 
ations on which the shorter-time-scale variations are later 
superimposed. For example, large-scale energy releases in 
the corona (e.g. through magnetic reconnection) might fur- 
ther sub-divide into a fractal structure, where the energy 
emitted by each sub-unit is proportional to the energy of 
its parent unit, but the time-scales for energy emission and 
the number of sub-units remain independent of total energy 
content (perhaps related to characteristic time-scales in the 
corona or disk). 

Interestingly, a linear dependence of RMS variability on flux 
is a natural outcome of the mechanism for producing red- 
noise variability proposed by Lyubarskii (1997), where fluc- 
tuations in mass accretion rate at different disk radii are 
propagated through the disk to produce variations on all 
time-scales in the inner disk (and associated corona). In this 
model, fractional variations in mass accretion rate are pro- 



5 CONCLUSIONS 

We have shown that the amplitude of the broadband noise 
variability in the lightcurves of the black hole candidate 
Cyg X-l and the millisecond X-ray pulsar SAX J1808.4-3658 
is dependent on flux, in that the RMS variability for a given 
segment of the lightcurve scales (on average) linearly with 
the segment mean flux. The linear relation between flux and 
RMS variability has a positive offset on the flux axis, sug- 
gesting the existence of a second, constant-flux component 
to the lightcurve which contributes ~ 25% of the total flux. 
The shape of the power spectrum (at least on short time- 
scales) remains independent of flux. The X-ray lightcurves 
of AGN are also consistent with the same linear scaling of 
RMS variability with flux, suggesting that this behaviour is 
an intrinsic feature of the broadband noise, which is char- 
acteristic of the lightcurves of compact accreting systems 
across at least 6 decades of central object mass. 
The spectrum of the constant component to the lightcurve 
is similar to the total spectrum, suggesting that it may cor- 
respond to quiet, non-variable regions in the X-ray emitting 
corona. The unusual behaviour of the variable component of 
the lightcurve, is inherently difficult to explain using stan- 
dard shot-noise models. An alternative, 'top down' approach 
is more suitable, where the longest-time-scale variations pre- 
cede the smallest. Possibilities include large coronal flaring 
regions which break down into a fractal structure of smaller 
sub-units, or models where variability is caused by variations 
in mass accretion rate which propagate through the disk, so 
that shorter time-scale variations from the inner parts of the 
disk are superimposed on longer time-scale variations from 
further out. 

If the linear RMS-flux relation is indeed common to the 
broadband noise variability of all compact accreting sys- 
tems, it should provide a useful tool for separating out com- 
ponents of lightcurves which are not related to the broad- 
band noise component. It will also prove interesting to see 
how other features of XRB power spectra (such as QPOs) 
scale with flux. 
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